Ergoregion instability of black hole mimickers* 
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Ultra-compact, horizonless objects such as gravastars, boson stars, wormholes and superspinars 
can mimick most of the properties of black holes. Here we show that these "black hole mimickers" 
will most likely develop a strong ergoregion instability when rapidly spinning. Instability timescales 
range between ~ 10 _5 s and ~ weeks depending on the object, its mass and its angular momentum. 
For a wide range of parameters the instability is truly effective. This provides a strong indication 
that astrophysical ultra-compact objects with large rotation are black holes. 



I. INTRODUCTION 



Black holes (BHs) in Einstein-Maxwell theory are char- 
acterized by three parameters [1]: mass M, electric 
charge Q and angular momentum J = aM ^ M 2 . 
BHs are thought to be abundant objects in the Uni- 
verse. Their mass is estimated to vary between 3Mq 
and 10 9 5 M Q or higher Q, their electrical charge is neg- 
ligible because of the effect of surrounding plasma [H and 
their angular momentum is expected to be close to the 
extremal limit because of accretion and merger events 
[J . A non-cornprehensive list of some astrophysical BH 
candidates 



comprener 

afia 



is shown in Table Q] 



Despite the wealth of circumstantial evidence, there 
is no definite observational proof of the existence of as- 
trophysical BHs due to the difficulty to detect an event 
horizon in astrophysical BH candidates [2|, |8| • Thus as- 
trophysical objects without event horizon, yet observa- 
tionally indistinguishable from BHs, cannot be excluded 
a priori. Some of the most viable alternative models 
describing an ultra-compact astrophysical object include 
gravastars, boson stars, wormholes and superspinars. 
Dark energy stars or gravastars are compact objects with 
de Sitter interior and Schwarzschild exterior Q. These 



two regions are glued together by a model-dependent in- 
termediate region. In the original model [9( the inter- 
mediate region is an ultra-stiff thin shell. Models with- 
out shells or discontinuities have also been investigated 

GOT. 

Boson stars are macroscopic quantum states which are 
prevented from undergoing complete gravitational col- 
lapse by Heisenberg uncertainty principle [l2j| . Their 
models differ in the scalar self-interaction potential which 
also set the allowed maximum compactness for a boson 
star. 

An exhaustive description of wormholes can be found in 
the monograph [u| (see also Ref . ) ■ In this work we 
shall consider particular wormholes which are infinitesi- 
mal variations of BH spacetimes. These wormholes may 
be indistinguishable from ordinary BHs (l5| . 
Superspinars are solutions of the gravitational field equa- 
tions that violate the Kerr bound. These geometries 
could be created by high energy corrections to Einstein 
gra vity such as those present in string-inspired models 



TABLE I: Mass, M, radius, R, angular momentum, J, and 
compactness, fi — M/R, for some BH candidates (from 0, [H, 
Mass and radius are in solar units. 
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Candidate 


M 


R x 


10~ b 


JjM 1 


/" = 


M/R 


GRO J 1655-40 


6.3 


1.6- 


- 2.6 


0.65 


-0.80 


0.47 


-0.83 


XTE J1550-564 


10 


2.1 - 


- 8.4 


0.90 


- 1.00 


0.25 


-0.99 


GRS 1915+105 


14 


2.9 - 


- 9.7 


0.98 


- 1.00 


0.30 


-0.99 


SGR A* 


4 x 10 f 


; < 


27 


0.50 


- 1.00 


> 


0.31 
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The objects described above can be almost as com- 
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pact as a BH and thus they are virtually indistinguish- 
able from BHs in the Newtonian regime, hence the name 
"BH mimickers" . Although exotic these objects provide 
viable alternatives to astrophysical BHs. BH mimick- 
ers being horizonless, no information loss paradox [17| 
arises in these spacetimes. Moreover they can be regu- 
lar at the origin, avoiding the problem of singularities. 
By Birkhoff 's theorem, the vacuum exterior of a spheri- 
cally symmetric object is described by the Schwarzschild 
spacetime. Thus the motion of orbiting objects both 
around a static BH and around a static ultra-compact 
object is the same and it makes virtually impossible to 
discern between a Schwarzschild BH and a static neutral 
BH mimickcr. Instead for rotating objects deviations in 
the properties of orbiting objects occur. Since BH mim- 
ickers are very compact these deviations occur close to 
the horizon and are not easily detectable elect romagnet- 
ically. To ascertain the true nature of ultra-compact ob- 
jects it is thus important to devise observational tests to 
distinguish rotating BH mimickers from ordinary Kerr 
BHs. The traditional way to distinguish a BH from a 
neutron star is to measure its mass. If the latter is larger 
than the Chandrasekhar limit, the object is believed to 
be a BH. However, this method cannot be used for the 
BH mimickers discussed above, because of their broad 
mass spectrum. The main difference between a BH and 
a BH mimicker is the presence of an event horizon in the 
former. Some indirect experimental methods to detect 
the event horizon has been proposed [HI, . Another 
very promising observational method to probe the struc- 
ture of ultra-compact objects is gravitational wave as- 
tronomy. From the gravitational waveform it is expected 
to detect the presence of an event horizon in the source 
[20| . Some other BH mimickers (for example electrically 
charged quasi-BHs [2l| ) are already ruled out by exper- 
iments. Moreover there are evidences that some model 
for BH mimickers is plagued by a singular behavior in 
the near- horizon limit [22|. 

Here, we describe a method originally proposed in 
[23l [2^ | for discriminating rotating BH mimickers from 
ordinary BHs. This method uses the fact that compact 
rotating objects without event horizon are unstable when 
an ergoregion is present. This ergoregion instability ap- 
pears in any system with ergoregions and no horizons 
[25| . The origin of this instability can be traced back to 
superradiant scattering. In a scattering process, super- 
radiance occurs when scattered waves have amplitudes 
larger than incident waves. This leads to extraction of 
energy from the scattering body [III [27], . Instability 
may arise whenever this process is allowed to repeat itself 
ad infinitum. This happens, for example, when a BH is 
surrounded by a "mirror" that scatters the superradiant 
wave back to the horizon, ampli fying it at each scatter- 
ing, as in the BH bomb process [2917 l30l] . If the mirror 
is inside the ergoregion, superradiance may lead to an 
inverted BH bomb. Some superradiant waves escape to 
infinity carrying positive energy, causing the energy in- 
side the ergoregion to decrease and eventually generating 



an instability. This may occur for any rotating star with 
an ergoregion: the mirror can be either its surface or, 
for a star made of matter non-interacting with the wave, 
its center. On the other hand BHs could be stable due 
to the absorption by the event horizon being larger than 
superradiant amplification. Indeed Kerr BHs are stable 
aganist small scalar, electromagnetic and gravitational 
perturbations [3l| . 

Rapidly rotating stars do possess an ergoregion and 
thus they are unstable. However typical instability 
timescales are shown to be larger than the Hubble time 
f32| . Thus the ergoregion instability is too weak to pro- 
duce any effect on the evolution of stars. This conclusion 
changes drastically for BH mimickers due to their com- 
pactness [111, [24| • For some of the rotating BH mimick- 
ers described above, instability timescales range between 
~ 10~ 5 s and ~ weeks depending on the object, its mass 
and its angular momentum. 

This paper is organized as follows. In Section [TT] we 
deal with gravastars and boson stars. We describe rotat- 
ing models for these objects and discuss their instability 
timescale. In Section IIIII a toy model for both rotat- 
ing wormholes and superspinars is presented. Section 
IIVI contains a brief discussion of the results and con- 
cludes the paper. Throughout the paper geometrizcd 
units (G = c = 1) are used, except during the discus- 
sion of results for rotating boson stars when we set the 
Newton constant to be G = 0.05/ (4tt) as in Ref. [33j |. 



II. GRAVASTARS AND BOSON STARS 

This section discusses the main properties of gravas- 
tars and boson stars as well as the method to compute 
the ergoregion instability for these objects. For a more 
detailed discussion see 1231. 



A. Nonrotating Gravastars 

Although exact solutions for spinning gravastars are 
not known, they can be studied in the limit of slow rota- 
tion by perturbing the nonrotating solutions [34j . This 
procedure was used in Ref. [35[ to study the existence 
of ergoregions for ordinary rotating stars with uniform 
density. In the following, we omit the discussion for the 
original thin-shell model by Mazur and Mottola @ and 
we focus on the anisotropic fluid model by Chirenti and 
Rezzolla QUEH. 

The model assumes a thick shell with continuous pro- 
file of anisotropic pressure to avoid the introduction of 
an infinitesimally thin shell. The stress-energy tensor is 
T^ v = diag[— p,p r ,pt,pt], where p r and pt are the ra- 
dial and tangential pressures, respectively. The spherical 
symmetric metric is 

ds 2 = -f(r)dt 2 + B(r)dr 2 + r 2 d^l (2.1) 
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and it consists of three regions: an interior (r < ri) 
described by a de Sitter metric, an exterior (r > r 2 ) 
described by the Schwarzschild metric and a model- 
dependent intermediate (ri < r < r 2 ) region. In the 
following we shall indicate with S — r 2 — T\ the thickness 
of the intermediate region and with p — M/r 2 the com- 
pactness of the gravastar. In the model by Chirenti and 
Rczzolla the density function is 

!po , < r < r\ interior 

ar 3 + br 2 + cr + d , n < r < r 2 intermediate 
, r 2 < r exterior 

where a, b, c and d are found imposing continuity condi- 
tions p(0) = p{ri) = p , p(r 2 ) = p'{r x ) = p'(r 2 ) = and 
po is found fixing the total mass, M. The metric coeffi- 
cients are 

'-I-***, (2.2) 
\ r 2 J B r 

where 

/ \ f r a 2 . Tit \ [ r 2m(r) + 87rr 3 p r 

m(r) = / 4?rr pdr , T(r) = / — , dr . 

Jo Jo r l r — ^m^rjj 

(2.3) 

The above equations and some closure relation, p r — 
p r (p), completely determine the structure of the gravas- 
tar [T(|. The behaviors of the metric coefficients for a 
typical gravastar are shown in Fig. [T] 



I 1 I 1 I 1 I 1 I 
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FIG. 1: Metric coefficients for the anisotropic pressure model 
(r a = 2.2, n = 1.8 and M = 1). 



of order 17, g t( p = —ujg^, where <p is the azimuthal co- 
ordinate and uj = u)(r) is the angular velocity of frame 
dragging. The full metric is 

ds 2 = -fdt 2 + Bdr 2 + r 2 d0 2 + r 2 sin 2 6 (d<j> - ujdtf , 

where /, B and u> are radial functions. If the gravastar 
rotates rigidly, i.e. f2 = constant, from the (t, <fi) compo- 
nent of Einstein equations we find a differential equation 
for w(r) H 

uj" + J + jj = 16nB(r) (w - ft) (p + p t ) , (2.5) 

where j = {fB) -1 / 2 is evaluated at zeroth order and p, 
Pt are given in terms of the nonrotating geometry. The 
above equation reduces to the corresponding equation 
for isotropic fluids (34j. Solutions of Eq. (|2.5D describe 
rotating gravastars to first order in SI. 

The ergoregion can be found by computing the surface 
on which g u vanishes [35[ ■ An approximated relation for 
the location of the ergoregion in very compact gravastars 
is 

= -fir) + Lu 2 r 2 sm 2 9. (2.6) 

The existence and the boundaries of the ergoregions can 
be computed from the above equations. We integrate 
equation (|2.5p from the origin with initial conditions 
(SI — uj)' — and (SI — uj) finite. The exterior solution 
satisfies lo = 2 J/r 3 , where J is the angular momentum of 
the gravastar. Demanding the continuity of both (CI — u))' 
and (A — u>), A and J are uniquely determined. The rota- 
tion parameter A depends on the initial condition at the 
origin. Figure [2] shows the results the gravastar model 




1. Slowly rotating gravastars and ergoregions 

Slowly rotating solutions can be obtained using the 
method developed in Ref. [3J|. A rotation of order A 
gives corrections of order A 2 in the diagonal coefficients 
of the metric (|2.ip and introduces a non-diagonal term 



FIG. 2: J I'M 1 and angular frequency fi for the anisotropic 
pressure model with r^ = 2.2, n = 1.8 and M = 1. 

described in the previous sections. The ergoregion can 
be located by drawing an horizontal line at the desired 
value of J/M 2 . The minimum of the curve is the mini- 
mum values of J/M 2 which are required for the existence 



4 



of the ergoregion. Comparison with the results for stars 
of uniform density [35j , shows that ergoregions form more 
easily around gravastars due to their higher compactness. 
The slow-rotation approximation is considered valid for 
VL/VLk < 1 where Mfljc = [i 3 / 2 is the Keplerian fre- 
quency. 

Depending on the compactness, ^i, the angular momen- 
tum, J, and the thickness, 5, a spinning gravastar does or 
does not develop an ergoregion. The formation of an er- 
goregion for rotating gravastar is exhaustively discussed 
in the whole parameters space in Ref. [36[. A delicate 
issue is the strong dependence on the thickness, 8, which 
cannot be directly measured by experiments. Figure [3] 
shows how the ergoregion width is sensitive to 8. 




FIG. 3: Ergoregion width (in units of M) as function of the 
thickness, 5 = r% — n, for T2 = 2.3, M = 1 and for different J 
values. From top to bottom: J/M 2 = 0.95, 0.90, 0.85, 0.80, 
0.75, 0.70, 0.65 and 0.60. The ergoregion width decreases as 
5 -v 0. 



symmetry otherwise. Since the Lagrangian density is in- 
variant under a global U(l) transformation, the current, 
= —i^d^Q+cc, is conserved and it is associated to a 
charge Q, satisfying the quantization condition with the 
angular momentum J = nQ [371 ] . The numerical proce- 
dure to extract the metric and the scalar field is described 
in Ref. [HJ ■ Throughout the paper we will consider solu- 
tions with n = 2, b = 1.1, A = 1.0, a = 2.0 and different 
values of (J,M) corresponding to J/{GM 2 ) ~ 0.566, 
0.731 and 0.858. In Fig. [H the metric functions for a bo- 
son star along the equatorial plane are shown. By com- 
puting the coefficient gu one can prove that boson stars 
develop ergoregions deeply inside the star. For this par- 
ticular choice of parameters, the ergoregion extends from 
r/{GM) — 0.0471 to 0.770. A more complete discussion 
on the ergoregions of rotating boson stars can be found 
in Ref. 




r/(GM) 



B. Rotating boson stars 

A example of rotating boson star is the model by Klei- 
haus, Kunz, List and Schaffer (KKLS) [Hj]. The KKLS 
solution is based on the Lagrangian for a self-interacting 
complex scalar field 

£kkls = -\<r + - vm) , (2.7) 

where U(\$\) = A|$| 2 (|$| 4 - a|$| 2 + 6). The mass of 
the boson is given by me = vA6. The ansatz for the 
axisymmetric spacetime is 




ds 2 = -fdr 



kg 

T 



dr 2 



d& 1 



r 2 sin 2 I 



9 



(dcp-ar)dty 



and $ = 6 e iu ° t+inv 



(2.8) 

where the metric components and 



the real function 4> depend only on r and 9. The require- 
ment that $ is single- valued implies n — 0, ±1, ±2, 

The solution has spherical symmetry for n = and axial 



FIG. 4: Left panel: Metric coefficients for a rotating bo- 
son star along the equatorial plane, with parameters n — 2, 



1.1, A = 1.0, a = 2.0, J/(GM 2 



0.566. Right panel: 



Fractional difference of the metric potentials between 9 = n/2 
and 6 — 7r/4 for the same star. 
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C. Ergoregion instability for rotating gravastars 
and boson stars 

The stability of gravastars and boson stars can be stud- 
ied perturbatively by considering small deviations around 
equilibrium. Due to the difficulty of handling gravita- 
tional perturbations for rotating objects, the calculations 
below are mostly restricted to scalar perturbations. How- 



ever the equation for axial gravitational perturbations of 
gravastars is identical to the equation for scalar pertur- 
bations in the large I = m limit [23| . There are also 
generic arguments suggesting that the timescale of grav- 
itational perturbations is smaller than the timescale of 
scalar perturbations for low m [38j |. Thus, scalar pertur- 
bations should provide a lower bound on the strength of 
the instability. 



TABLE II: WKB results for the instability of rotating gravastars with V2 = 2.2, ri = 1.8 and M = 1. 



t/M 


in 


J/M 1 = 0.40 JIM 1 = 0.60 J/M' 2 = 0.80 J/M 2 = 0.90 J/M 2 = 1.0 
Q/Q K = 0.33 Q/Q K = 0.49 tl/{l K = 0.65 Q/Q K = 0.74 Q,/£l K = 0.82 


1 


1.33 x 10 7 


2.78 x 10 4 


5.99 x 10 3 


3.58 x 10 3 


2.34 x 10 3 


2 


8.25 x 10 7 


1.14 x 10 6 


1.11 x 10 5 


4.81 x 10 4 


2.33 x 10 4 


3 


1.31 x 10 10 


5.65 x 10 7 


2.25 x 10 6 


6.82 x 10 5 


2.45 x 10 5 


4 


2.50 x 10 12 


2.95 x 10 9 


4.81 x 10 7 


1.02 x 10 7 


2.73 x 10 6 


5 


5.06 x 10 14 


1.59 x 10 11 


1.02 x 10 9 


1.52 x 10 8 


3.07 x 10 7 



1. Scalar field instability for slowly rotating gravastars: 
WKB approach 

Consider now a minimally coupled scalar field in the 
background of a gravastar. The metric of gravastars is 
given by Eq. (|2.4p . In the large / = m limit, which is 
appropriate for a WKB analysis [H, H^] , the scalar field 
can be expanded in spherical armonics, Yi m — Yi m (9 , <f) 
as 



$ = ^2xim exp 

lm 



dr 



-iu)t-y 



2/ 2B / 

(2.9) 

The functions xim = Xim( r ) are determined by the 
Klein-Gordon equation which, dropping terms of order 
O (1/m 2 ), yields 



Xim + m 2 T{r , S)x(m = , 



(2.10) 



where E = —uj/m and 



T 



B(r) 
f(r) 



(E - V+) (E - VL) , V± 



Equation (f2~T0l) can be shown to be identical for the axial 
gravitational perturbations of perfect fluid stars [24[ • 

The WKB method [13] for computing the eigenfre- 
quencies of Eq. (|2.10p is in excellent agreement with full 
numerical results [3!|. The quasi-bound unstable modes 
are determined by 



/?? 



0,1,2,. 



(2.11) 



and have an instability timescale 



4 exp 



2m 



T\dr 



d 



Tdr. 



(2.12) 



where r a , are solutions of V+ = E and r c is determined 
by the condition V- — E. 

Table [III shows the WKB results for the anisotropic 
pressure model for different values of J/M 2 . Although 
the WKB approximation breaks down at low m values, 
these results still provide reliable estimates [Hf. This 
claim has be verified with a full numerical integration of 
the Klein-Gordon equation. The results show that the 
instability timescale decreases as the star becomes more 
compact. Larger values of J/M 2 make the star more un- 
stable. The maximum growth time of the instability can 
be of the order of a few thousand M, but it crucially de- 
pends on J, fi and S 36] . For a large range of parameters 
this instability is crucial for the star evolution. Gravita- 
tional perturbations are expected to be more unstable. 
Moreover it is worth to notice that the slowly rotating 
approximation allows only for /i < 0.5, while for rotating 
BHs 0.5 < fi < 1 (see Tabic fl}. The ergoregion insta- 
bility being monotonically increasing with /i, we expect 
that instability timescales for realistic gravastars should 
be much shorter than the ones computed. For most of the 
BH mimickers models to be viable we require J/M 2 ~ 1 
and fj, ~ 1. ft would be interesting to study whether 
the ergoregion instability is or is not always effective in 
this case. Possible future developments include: (i) a full 
rotating gravastar model, which allows for /i > 0.5; (ii) 
the stability analysis against gravitational perturbations 
for rotating gravastars; (iii) a gravavastar model which is 
not strongly dependent on the thickness, d. 
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The ergoregion instability of a rotating boson star 
is straightforwardly computed following the method de- 
scribed above for spinning gravastars. We refer the reader 
to [23| and we only summarize the results in Table IIIII 
The maximum growth time for this boson star model is 
of the order of 10 6 Af for J/GM 2 = 0.857658. Thus the 
instability seems to be truly effective for rotating boson 
stars. 



TABLE III: Instability for rotating boson stars with parame- 
ters n — 2, b = 1.1, A = 1.0, a = 2.0 and different values of J 
(from [231]). The Newton constant is denned as 4nG = 0.05. 



t/(GM) 


m 


J/GhP = 0.5661 J /GAP = 0.7307 J /GAP = 0.8577 


1 


8.847 x 10 2 


6.303 x 10 3 




2 


7.057 x 10 3 


5.839 x 10 4 


1.478 x 10 6 


3 


6.274 x 10 4 


9.274. x 10 s 


2.815 x 10 8 


4 


5.824 x 10 5 


1.603 x 10 7 


2.815 x 10 10 


5 


5.554 x 10 6 


2.915 x 10 8 


1.717 x 10 12 



A rigorous analysis of the ergoregion instability for these 
models is a non-trivial task. Indeed known wormhole 
solutions are special non-vacuum solutions of the gravi- 
tational field equations, thus their investigation requires 
a case-by-case analysis of the stress-energy tensor. More- 
over exact solutions of four-dimensional superspinars are 
not known. To overcome these difficulties, the follow- 
ing analysis will focus on a simple model which captures 
the essential features of most Kerr-like horizonless ultra- 
compact objects. Superspinars and rotating wormholes 
will be modeled by the exterior Kerr metric down to their 
surface, where mirror-like boundary conditions are im- 
posed. This problem is very similar to Press and Teukol- 
sky's "BH bomb" [H, [3(3], i.e. a rotating BH surrounded 
by a perfectly reflecting mirror with its horizon replaced 
by a reflecting surface. For a more detailed discussion 
see 
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III. A TOY MODEL FOR KERR-LIKE OBJECTS 



2. Superspinars and Kerr-like wormholes 



This section discusses Kerr-like objects such as partic- 
ular solutions of rotating wormholes and superspinars. 



A superspinar of mass M and angular momentum J 
aM can be modeled by the Kerr geometry [161 ] 



J 



ds Kc 



2Mr 



dt 2 + ^-dr 2 
A 



(r 2 + a 2 ) 2Mr 



sm 



+ 



sm 



4Mr 



a sin 2 Qd(j>dt + £d6> 2 



r 



(3.1) 



where £ = r 2 + a 2 cos 2 8 and A = r 2 + a 2 — 2Mr. Un- 
like Kerr BHs, superspinars have a > M and no horizon. 
Since the domain of interest is — oo < r < +oo, the 
space-time possesses naked singularities and closed time- 
like curves in regions where g^ < (4p|. High energy 
modifications (i.e. stringy corrections) in the vicinity of 
the singularity are also expected. 

Kerr-like wormholes are described by the metric 



ds 



wormhole 



dsT, 



5g a bdx a dx b 



(3.2) 



where 8g a b is infinitesimal. In general, Eq. (13. 2j) de- 
scribes an horizonless object with a excision at some 
small distance of order e from the would-be horizon 
[i~5| . Wormholes require exotic matter and/or divergent 
stress tensors, thus some ultra-stiff matter is assumed 
close to the would-be horizon. In the following, both 
superspinars and wormholes will be modeled by the 



Kerr metric with a rigid "wall" at finite Boyer-Lindquist 
radius Vq, which excludes the pathological region. 

A. Instability analysis 



If the background geometry of superspinars and worm- 
holes is sufficiently close to the Kerr geometry, its per- 
turbations is determined by the equations of perturbed 
Kerr BHs [24{ . Thus the instability of superspinars and 
wormholes is studied by considering Kerr geometries with 
arbitrary rotation parameter a and a "mirror" at some 
Boyer-Lindquist radius rg. Using the Kinnersley tetrad 
and Boyer-Lindquist coordinates, it is possible to sepa- 
rate the angular variables from the radial ones, decou- 
pling all quantities. Small perturbations of a spin-s field 
are reduced to the radial and angular master equations 

m 



j 



A _ s d_ 
dr 



idR. 



I in 



dr 



K 2 - 2is(r - M)K 



+ Aisuor — A 



Rim — , 



(3.3) 
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[(1 - X 2 ) 3 Si m ,x] x + 



(auix) — 2au>sx + s 
I 



A-l m — 



(m - 



sx) 



sSlm — , 



(3.4) 



where x = cos 9, A = r 2 — 2Mr + a 2 and K = 
(r 2 + a 2 )uj — am. Scalar, electromagnetic and gravita- 
tional perturbations correspond to s — 0, ±1, ±2 respec- 
tively. The separation constants A and s Ai m are related 
by A = s Ai m + a 2 uj 2 — 2amuj. 



1. Analytic results 
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FIG. 5: Imaginary and real parts of the characteristic gravi- 
tational frequencies for an object with a = 0.998A/, according 
to the analytic calculation for rapidly-spinning objects. The 
mirror location is at ro = (1 + e)r+. The real part is ap- 
proximately constant and close to mQ, in agreement with the 
assumptions used in the analytic approach. 

Following Starobinsky [13], equations (|3.3p - (|3.4[) can 
be analytically solved in the slowly-rotating and low- 
frequency regime, loM <C 1, and in the rapidly-spinning 
regime, where r + ~ r_ and u> ~ mflh, where tth = 
a/(2Mr + ) is the angular velocity at the horizon. The 
details of the analytic approximation are described in 
Ref. (24|. Analytic solutions for a star with a = 0.998A/ 
are shown in Fig. [5] where gravitational perturbations are 
considered. The instability timescale for gravitational 
perturbations is about five orders of magnitude smaller 
than the instability timescale for scalar perturbations. 



B. Instability analysis: numerical results 

The oscillation frequencies of the modes can be found 
from the canonical form of Eq. 



where 



Y = A s / 2 (r 



d 2 Y 
dr 2 



) 1/2 R, 



VY = 0, 



(3.5) 



V 



K 2 - 2is(r - M)K + A(4irujs - A) 
(r 2 + a 2 ) 2 



G 



dG 

dr* 



and K = (r 2 + a 2 )uj - am, G = s(r - M)/(r 2 + a 2 ) + 



rA(r 2 



The separation constant A is related to 



the eigenvalues of the angular equation by A = s Ai m + 
2amio. The eigenvalues s Ai m are expanded in 



a 2 io 2 



power series of aw as [42| 

sA lm = J2f^) h 



(3.6) 



fc=0 



Terms up to order (auj) 2 are included in the calculation. 
Absence of ingoing waves at infinity implies 



(3.7) 



Y 



Numerical results are obtained by integrating Eq. (J375J) 
inward from a large distance r^. The integration is 
performed with the Runge-Kutta method with fixed cj 
starting at Mr oo = 400, where the asymptotic behav- 
ior (|3.7p is imposed. (Choosing a different initial point 
does not affect the final results.) The numerical integra- 
tion is stopped at the radius of the mirror ro, where the 
value of the field Y(oj,ro) is extracted. The integration 
is repeated for different values of u> until Y(ui,ro) = is 
obtained with the desired precision. If Y (lo, ro) vanishes, 
the field satisfies the boundary condition for perfect re- 
flection and to = luq is the oscillation frequency of the 
mode. 



Objects with a < M 



The regime a < M requires a surface or mirror at 
ro = + e) > r + . Thus the compactness is M/r$ ~ 
(1 — e)M/r + and, in the limit e — > 0, it is infinitesi- 
mally close to the compactness of a Kerr BH. Numerical 
results for scalar and gravitational perturbations of ob- 
jects with a < M are summarized in Table IIVI and are 
in agreement with the analytic results [24|. The insta- 
bility is weaker for larger m. This result holds also for 
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TABLE IV: Characteristic frequencies and instability 
timescales for a Kerr-like object with a = 0.998M. The mir- 
ror is located at e = 0.1, corresponding to the compactness 

fJ, ~ 0.9jLtKerr- 

(Re(tu)M ,lm(w)M) 
I — m s = s = 2 

1 (0.1120,0.6244 x 10" 5 ) 

2 (0.4440,0.5373 x 10" 5 ) (0.4342,0.2900) 

3 (0.7902,0.1928 x 10 -5 ) (0.7803,0.2977) 

4 (1.1436,0.5927 x 10~ 6 ) (1.1336,0.3035) 



/ 7^ m and s = 0, ±1 and ±2. The minimum instability 
timescale is of order r ~ 10 5 M for a wide range of mir- 
ror locations. Figure [6] shows the results for gravitational 
perturbations. Instability timescales are of the order of 
t ~ 2 6M. Thus gravitational perturbations lead to an 
instability about five orders of magnitude stronger than 
the instability due to scalar perturbations (see Table Hv|) . 
Figure [5] shows that the ergoregion instability remains 
relevant even for values of the angular momentum as low 
as a = Q.6M. 




FIG. 6: Details of the instability for gravitational perturbations, for different I = m modes and a/M = 0.998 (top panels) and 
for I — m — 2 and different a/M < 1. 



2. Objects with a> M 

Objects with a > M could potentially describe su- 
perspinars. Several arguments suggest that objects ro- 
tating above the Kerr bound are unstable. Firstly, ex- 
tremal Kerr BHs are marginally stable. Thus the ad- 
dition of extra rotation should lead to instability. Sec- 
ondly, fast-spinning objects usually take a pancake-like 
form [43j and are subject to the Gregory-Laflammc in- 
stability [44, 45]. Finally, Kerr-like geometries, like naked 
singularities, seem to be unstable against a certain class 
of gravitational perturbations [III [43] called algebraically 



special perturbations For objects with a > M the 

surface or mirror can be placed anywhere outside r = 0. 
In general the instability is as strong as in the a < M 
regime. An example in shown in Fig. [7] for the surface at 
tq/M = 0.001. This result confirms other investigations 
suggesting that ultra-compact objects rotating above the 
Kerr bound are unstable 14811. 
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FIG. 7: The fundamental I = m — 2, 3, 4 modes of an object 
spinning above the Kerr bound as function of rotation. The 
surface is located at ro/M = 0.001. 

IV. CONCLUSION 

We investigated the ergoregion instability of some 
ultra-compact, horizonless objects which can mimick the 
spacetime of a rotating black hole. We studied some of 
the most viable BH mimickers: gravastars, boson stars, 
wormhole and superspinars. 

If rotating, boson stars and gravastars may develop er- 
goregion instabilities. Analytical and numerical results 
indicate that these objects are unstable against scalar 
field perturbations for a large range of the parameters. 



Slowly rotating gravastars can develop an ergoregion de- 
pending on their angular momentum, their compactness 
and the thickness of their intermediate region. In a recent 
work [36| it has pointed out that slowly rotating gravas- 
tars may not develop an ergoregion. In the formation of 
the ergoregion for rotating gravastars an important role 
is played by the thickness (see Figure [3]) which is not eas- 
ily detectable. Thus further investigations are needed to 
better understand the ergoregion formation in physical 
resonable gravastar models. 

The instability timescale for both boson stars and 
gravastars can be many orders of magnitude stronger 
than the instability timescale for ordinary stars with uni- 
form density. In the large / = m approximation, suitable 
for a WKB treatment, gravitational and scalar pertur- 
bations have similar instability timescales. In the low-rn 
regime gravitational perturbations are expected to have 
even shorter instability timescales than scalar perturba- 
tions. Instability timescales can be as low as ~ 0.1 sec- 
onds for a M = IMq objects and about a week for su- 
permassive BHs, M — 10 6 M Q , monotonically decreasing 
for larger rotations and a larger compactness. 

The essential features of wormholes and superspinars 
have been captured by a simple model whose physical 
properties are largely independent from the dynamical 
details of the gravitational system. Numerical and ana- 
lytic results show that the ergoregion instability of these 
objects is extremely strong for any value of their angular 
momentum, with timescales of order 10~ 5 seconds for a 
IMq object and 10 seconds for a M — 10 6 M Q object. 
Therefore, high rotation is an indirect evidence for hori- 
zons. 

Although further studies are needed, the above investi- 
gation suggests that exotic objects without event horizon 
are likely to be ruled out as viable candidates for as- 
trophysical ultra-compact objects. This strengthens the 
role of BHs as candidates for astrophysical observations 
of rapidly spinning compact objects. 
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